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CHAPTER  1. 
INTRODUCTION 

1-1  General  Discussion 

A  vast  quantity  of  organic  compounds  has  been  introduced  recently  to 
improve  our  living  standards.   Most  of  these  compounds  are  the  products  from 
the  procedure  of  synthesis  by  industries.   The  discharge  of  synthetic  organic 
compounds  into  the  environment  will  cause  severe  adverse  effects  including 
toxicity,  carcinogencity,  taste  and  odor  problems,  and  degradation  of  the 
quality  of  water  for  consumptive  use.   The  potential  impacts  of  hazardous 
organic  pollutants  in  wastewater  constitute  a  matter  of  steadily  expanding 
concern  for  the  water  quality  specialist.   Control  of  toxic  pollutants  is  also 
gaining  increased  emphasis  in  both  water  and  wastewater  treatment.   Therefore. 
it  is  necessary  to  regulate  the  release  of  these  organics  to  the  environment 
and  find  methods  to  reduce  or  effectively  remove  these  contaminents  from 
water . 

The  Environmental  Protection  Agency  <EPA)  has  developed  a  list  of  harmful 
chemicals.   In  the  list,  129  pollutants  are  termed  as  priority  pollutants. 
The  reasons  for  singling  out  these  priority  pollutants  are  because  of  their  : 
(1)  frequencies  of  occurrence  in  water,  (2)  chemical  stability ,  (3)  quantities 
of  the  chemical  produced,  (4)  availability  of  chemical  analysis  methods. 

Phenol  and  its  derivatives  occupy  a  prominent  position  on  the  EPA  priority 
pollutant  list.   In  addition  to  being  suspected  carcinogens,  phenol  and  some 
of  its  derivatives  have  been  shown  to  be  either  toxic  or  lethal  to  fish  at 
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concentrations  of  5  to  25  mg/1  (Khan,  1981).   The  presence  of  phenol  at  a 
concentration  of  only  2  ug/1  imparts  objectionable  tastes  and  odors  to  drin- 
king water  when  combined  with  chlorite,  due  to  the  formation  of  chlorophenols. 
In  referring  to  the  phenolic  compounds,  it  is  intended  to  include  not  only 
phenol  itself  but  also  chlorinated  phenols.   As  a  concequence,  though  p- 
chlorophenol  (PCP)  is  not  found  in  the  EPA  list  of  priority  pollutants,  it  is 
of  concern  in  this  study  due  to  the  fact  that  it  is  a  member  of  the  chloro- 
phenol  family  which  was  included  in  the  list. 

Phenol  is  recovered  from  coal  tar.  and  large  amounts  are  produced  synthe- 
tically.  It  is  used  extensively  in  the  synthesis  of  organic  products,  parti- 
cularly in  the  production  of  phenolic  type  resins.   Sources  of  phenolic  com- 
pounds occurring  as  a  natural  compound  in  industrial  wastes  include  effluents 
from  coal-gas,  coal-coking,  and  petroleum  industries  as  well  as  in  a  wide 
variety  of  industrial  wastes  from  processes  involving  the  use  of  phenol  as  a 
raw  material.   Industries  having  the  chlorophenol  family  in  their  wastewater 
are  textile,  ink  metal  finishing,  steam  electric,  and  leather  tanning  and 
finishing,  etc. 

Removal  or  reduction  of  these  organic  compounds  to  an  acceptable  level  has 
become  an  increasing  concern  in  the  fields  of  water  and  wastewater  treatment 
technology.   Considerable  efforts  have  been  made  by  the  EPA  and  other  resear- 
chers to  improve  conventional  techniques.   The  occurrence  of  maximum  concen- 
trations of  phenol  and  p-chlorophenol  in  the  raw  wastewater  and  in  treated 
industrial  effluents  are  listed  in  Tables  1.1  and  1.2. 

In  the  realm  of  dissolved  organic  contaminent  control,  treatment  methodo- 
logies such  as  biological  degradation,  adsorption,  ion  exchange,  chemical 


oxidation,  membrane  separations,  and  incineration,  have  been  applied.   Adsor- 
ption on  activated  carbon  has  gained  wide  acceptance  and  is  regarded  as  the 
most  efficient  and  economical  procedure  for  recovering  undesirable  materials 
from  dilute  aqueous  solution.   Activated  carbon  is  obtained  through  a  control- 
led oxidation  process  which  results  in  a  porous  carbon  structure  with  a  large 
surface  area.   The  large  surface  area  gives  the  activated  carbon  a  high  capa- 
city to  adsorb  organic  materials. 

Adsorption  processes,  especially  those  that  use  granular  activated  carbon, 
will  find  increasing  use  in  wastewater  and  potable  water  treatment.   It  is 
necessary  to  know  the  adsorption  capacity  as  well  as  the  adsorption  rate  to 
design  adsorption  equipment.   Adsorption  capacity  is  usually  expressed  by  an 
isotherm  based  on  measured  data.   The  rate  of  adsorption  depends  on  both 
transport  rate  from  the  outer  surface  of  the  particles  and  also  on  the  tran- 
sport rate  into  the  particles.   One  or  both  of  these  rates  may  determine  the 
total  rate  of  adsorption. 

A  breakthrough  curve  from  the  operation  of  carbon  columns  can  be  used  to 
obtain  the  performance  of  activated  carbon  columns  to  remove  specific  organic 
pollutants  and  to  determine  the  contact  time  required  to  produce  the  desired 
effluent  concentration.  Since  operation  of  such  columns  is  time  consuming  and 
expensive,  several  predictive  mathematical  models  have  been  proposed  to  elimi- 
nate or  greatly  reduce  the  need  for  pilot  studies.   These  models  which  are 
used  for  forecasting  adsorption  column  require  reaction  coefficients  for  adsor- 
ption capacity  and  the  adsorption  rates.  The  former  is  obtained  from  batch 
equilibrium  system;  and  the  latter  is  obtained  from  batch  kinetic  system. 
These  coefficients  are  determined  from  the  aid  of  mathematical  models. 


Many  researchers  have  been  working  on  the  phenomenon  of  activated  carbon 
adsorption  and  have  proposed  many  mathematical  models  in  the  literature. 
Almost  all  of  them  are  based  on  the  constant  temperature  assumption,  while 
very  a'  few  investigations  incorporate  the  study  of  temperature  effect  and  then 
only  on  a  limited  amount  of  data.  It  is  noted  that  temperature  might  be 
changed  abruptly  or  unexpectedly  during  operation  under  pratical  conditions. 
The  main  purpose  of  this  research  is  to  take  into  account  the  effects  of 
variation  of  temperature  on  adsorption.   The  batch  kinetic  and  fixed  bed 
(column)  studies  were  conducted  with  constant  temperatures  and  compared  with 
step  change  in  temperature. 

It  will  be  helpful  if  a  model  incorporates  temperature  as  a  variable.   A 
mathematical  model  developed  to  compare  the  performance  of  a  fixed  bed  adsor- 
ber at  constant  temperture  was  modified  to  take  into  account  influent  temper- 
taure  variations.   Experiments  were  conducted  to  verify  the  predicted  values 
from  the  modified  mathematical  model. 


Table   1.1 

Reported  Maximum  Concentrations  for  Phenol 
in  Industrial  Effluents 


Industrial 
category 


Max .  Concentration 
Raw  water  <;*g/l) 


Reported 
Treated  (^jg/1) 


Paint  Formulating  3,800 

Gum  and  Wood  Chenicals  23,000 

Leather  Tanning  2,500 

Ink  Formulating  536 

Textile  Industry  4,900 

Petroleum  Refineing 

Pulp,  Paper  and  Paperboard 

Industry  640 

Coal  Mining  3 

Steam  Electric  Power 

Plant  100 

Pulp  and  Paperboard  Mills 

and  Paper  Coverts  624 

Organic  Chemicals 

Manufacturing  7,300 

Paramaceuticals  51,000 


1.240 
1.900 
1.400 

18 

50 

4.200 

250 


89 

30 
120 


Source  :  EPA  Development  Document  for  Effluent  Limitation  Series 


Table  1.2 

Reported  Maximum  Concentrations  for  Chlorophenol 
in  Industrial  Effluents 


Industrial 
category 


Max.  Concentration 

Reported 

Raw  water 

ilig/1) 

Treated  (/*g/l) 

4.900 

3,700 

47 

-- 

306 

-- 

1,300 

-- 

940 

15 

Paint  Formulating 

Gun  and  Wood  Chemcals 

Timber  Product 
Processing 

Ink  Formulating 

Textile  Hills 

Petroleum  Ref ineing 

Pulp,  Paper  and  Paperboard 
Industry 

Lether  Tanning  and 
Finishing 

Steam  Electric  Power 
Plant 

Pulp  and  Paperboard  Mills 
and  Paper  Coverts 

Organic  Chemicals 
Manufacturing 

Paramaceuticals 


1,200 


6,200 


12 


4,200 


1,400 


1,700 


62 


Source  :  EPA  Development  Document  for  Effluent  Limitation  Series 


1-2  Objectives  and  Scope  of  the  Research 
The  objectives  of  the  research  are: 

1.  To  establish  adsorption  isotherm  for  phenol  and  PCP  at  constant 
temperature  and  extend  an  available  isotherm  equation  for  various  tempera- 
tures. 

2.  To  assess  the  performance  of  an  adsorption  model  to  predict  adsorption 
rates  at  constant  temperture  and  step  change  in  temperature  from  single 
and  multicomponent  systems. 

3.  To  assess  the  ability  of  a  fixed  bed  model  to  predict  breakthrough  curves 
at  constant  and  step  change  in  temperature  from  single  component  system. 

4.  To  compare  the  breakthrough  curves  for  step  changes  in  temperaure  for 
multicomponent  systems. 

The  scopes  of  this  research  are  as  follow: 

1.  Establishment  of  adsorption  equilibrium  and  kinetics  of  phenol  and  PCP. 

(1)  Single  solute  and  multicomponent  equilibrium  studies  were  conducted  to 
obtain  isotherms  at  four  different  temperatures. 

(2)  Batch  kinetic  studies  were  conducted  to  determine  adsoprtion  rates  and 
transport  parameters .   Adsorption  rates  were  developed  for  single 
solute  and  multicomponent  systems  for  constant  temperature  and  for 
step  changes  in  temperature. 

2.  Fixed  bed  (column)  studies  were  conducted  at  two  different  temperatures . 
The  experimental  data  from  single  component  tests  were  used  to  compare 


with  prediction  using  a  mathematical  model.   The  goal  was  to  show  the 
validity  of  the  model  in  predicting  the  breakthrough  curves  of  fixed  bed 
adsorption.   For  multicomponent  tests,  the  results  were  applied  to  compare 
adsorption  ability  with  different  procedures  of  step  change  in  tempera- 
tures. 


CHAPTER   2. 
REVIEW  OF  RELATED  LITERATURE 

2-1   Single  Solute  Adsorption  Models 

2-1-1   Equilibrium  Models 

Data  collected  during  an  equlibrium  test  will  describe  the  performance  of 
the  activated  carbon.   Several  mathematical  relationships  have  been  developed 
to  describe  the  equilibrium  distribution  of  solute  between  the  solid  and 
liquid  phases.   The  equilibrium  relationships  defined  as  adsorption  isotherm 
are  plots  between  equilibrium  concentration.  Ce,  which  is  the  amount  of  orga- 
nic compound  left  in  solution,  and  surface  concentration,  qe,  the  amount  of 
compound  on  the  surface  of  the  activated  carbon. 

An  adsorption  isotherm  equation  usually  represents  equilibrium  data  at  a 

given  temperature.   Three  of  the  most  common  isotherm  equations  are  the 

Langmuir  isotherm,  the  Freundlich  isotherm,  and  the  Brunauer-Emmett-Teller 

(BET)  isotherm.   The  Freundlich  isotherm  (Freundlich,  1926)  which  is  basically 

empirical,  is  represented  as 

qe  =  K  0,1'"  (2.1) 

or 

log  qe  ■  log  K  *  1/n  log  Ce  (2.2) 

A  plot  of  log  qa  vs.  log  CQ  should  yield  a  straight  line  for  adsorption  data 
which  follow  the  Freundlich  theory.   The  values  of  the  constants  K  and  n  that 
must  be  evaluated  for  each  solute  and  temperature  are  easily  obtained  from  the 
straight  line.   Although  the  assumption  was  based  on  the  fact  that  the  adsor- 
bent had  a  heterogeneous  surface  composed  of  different  classes  of  adsorption 
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sites,  ths  Freundlich  isotherm  is  often  useful  as  a  leans  for  data  descrip- 
tion.  The  Langauir  isothera  (Langmuir,  1318)  is  based  on  the  assumption  that 
the  adsorbed  layer  will  be  aonoaolecular.   The  Langauir  isothera  is  coaaonly 
written  as 

0  b  C, 
qa  •   (2.3) 

1'bC, 

In  addition  to  the  previous  assumption,  the  Langauir  isothem  also  assumes 
that  energy  of  adsorption  is  uniform.  Thus,  0  is  the  number  of  soles  of  solute 
adsorbed  by  per  unit  mass  of  adsorbent  in  forming  a  complete  monolayer  on  the 
adsorbent.   The  constant,  b,  is  related  to  the  energy  of  enthalpy  of  adsor- 
ption.  Taking  the  reciprocal  of  both  sides  of  Equation  2.3  yields  two  conve- 
nient linear  forms  of  the  Langauir  equation. 
C.         1  Ce 


qe        Ob  0 

1  1  1 


(2.4) 


(2.5) 


qe         0  0  b    Ce 

A  straight  line  will  be  obtained  when  the  quantity  1/qe  is  plotted  against 
1/Ce.   Values  of  the  constant,  Q  and  b,  can  then  be  determined  from  the  slope 
and  intercept  of  the  plot. 

Neither  the  Freundlich  nor  Languair  isothera  aay  describe  the  data  satis- 
factorily over  a  wide  range  of  concentration;  therefore,  an  eapirical  equation 
with  three  parameters  was  applied  by  Redlich  and  Peterson  (1959),  Radke  and 
Prausnitz  (1972),  and  Mathews  and  Weber  (1977)  to  describe  adsorption  isotherm 
when  the  concentration  of  the  adsorbent  varied  widely. 

A   Ce 

qa  •   fl  <  1  (2.6) 

1  *  B  0/ 
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At  low  concentrations,  this  three  parameter  equation  reduced  to  a  linear  form: 
while  at  high  concentration,  it  becomes  the  Freundlich  isotherm.   For  the 
special  case  of  $   -    1  ,  it  becomes  the  Langmuir  isotherm.   The  parameters,  A, 
B  ,  and  &   have  to  be  determined  by  the  best  statistical  representation  of  the 
experimental  data.   In  this  research,  the  three  parameter  equation  is  selected 
to  describe  equilibrium  adsorption. 

BET  isotherm  based  on  the  assumption  that  molecules  could  be  adsorbed  more 
than  one  layer  thick  on  the  surface  of  the  adsorbent.  The  equation  is  commonly 
written  as 

qe  =  (2.7) 

<CS-C>  C  1  ♦  (A-l)  C/Cs  3 

where  qe  =  amount  of  the  solute  adsorbed  per  unit  mass  of  adsorbent 

C  =  concentration  of  the  solute  in  solution  at  equilibrium 

C6  =  saturation  concentration  of  the  solute 

XB  =  amount  of  the  solute  adsorbed  in  forming  a  complete  monolayer  ■ 

A   -  a  constant  to  describe  the  energy  of  interaction  between  the 

the  solute  and  the  adsorbent  surface 
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2-1-2  Batch  Kinetic  Models 

The  equilibrium  distribution  of  solute  between  the  liquid  and  solid  phases 
is  an  important  property  of  an  adsorption  system.   Of  equal  importance  to  the 
engineer  is  the  kinetics  of  the  system  which  describe  the  rate  at  which 
equilibrium  is  reached.   Process  kinetics  describing  the  rates  at  which  mole- 
cules are  transferred  from  solution  to  the  surface  of  the  carbon  particles 
involve  several  sequential  and  parallel  transport  and  reaction  phenomena. 
Four  distinct  steps  must  take  place  for  adsorption  to  occur: 

(1)  The  adsorption  solute  passes  through  a  film  surrounding  the  adsorbent 
particle  to  the  surface  of  the  particle,  (film  transfer) 

(2)  The  adsorbed  solute  be  transferred  to  an  adsorption  site  on  the  inside 
of  pore,  (pore  diffusion) 

<3)  Ths  adsorbed  solute  diffuses  along  the  surface  of  the  pore,  (surface 

diffusion) 
(4)  The  solute  becomes  attached  to  the  surface  of  the  pore,  i.e.,  be 
adsorbed .  ( adsorption  > 
Figure  2.1  represents  the  overall  rates  of  adsorption  of  solute  onto  granular 
particles.   Many  factors  influence  the  rates  at  which  adsorption  occur  and  the 
extent  to  which  a  particular  material  can  be  adsorbed.   Several  of  the  factors 
will  be  discussed  in  later  sections. 

There  are  many  models  which  are  commonly  used  for  describing  one  or  a 
combination  of  different  rate  processes.   Differnent  assumptions  were  made  to 
simplify  rate  process  by  these  models.   Four  of  these  models  are  : 

A  second-order  reaction  rate  model  was  developed  by  Thomas  (1949)  for 
adsorption  and  ion-exchange  reactions.   Hiester  and  Vermuelen  (1952)  extended 
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it  to  diffusion  controlled  adsorption.   Keinath  and  WeDer  used  this  method  to 
predict  breakthrough  curves  for  expanded  fixed  bed  adsorption.   In  this 
method,  the  adsorption  rate  is  equated  to  the  diffusional  rates  to  obtain  the 
reaction  rate  constant  in  terns  of  the  transport  rate  constant. 

A  second  type  of  modeling  approach  is  the  apDlication  of  linear  driving 
force  expressions  for  internal  and  external  transport  processes  (Glueckauf . 
1955;  Wilde,  I960;  Tien.  1981). 

The  pore  phase  diffusion  model  considers  the  adsorbent  particles  as 
consisting  of  a  solid  phase  interspersed  with  very  small  pores.   The  adsorbate 
diffuses  into  the  pores  in  the  fluid  phase,  and  adsorption  occurs  at  the 
internal  surface  (Edeskuty  and  Amundson.  1952;  Kasten  .  1952:  DiGiano  and 
Weber.  1969,  1972). 

The  homogeneous  solid  phase  diffusion  model  proposed  by  Rosen  (1952)  will 
be  outlined  in  detail.   The  model  is  developed  for  predicting  adsorption  rates 
in  batch  and  fixed  bed  reactors .  The  model  equations  are  solved  numerically 
for  three  parameter  isotherm  with  the  film  transfer  and  solid  phase  diffusion 
as  the  rate  controlling  steps.   In  a  rapidly  agitated  batch  reactor,  the 
concentration  of  the  solute  C  and  the  distribution  of  the  adsorbent  particle 
may  be  assumed  to  be  uniform  throughout  the  reactor.   The  film  transfer  from 
the  fluid  phase  to  the  solid  surface  is  expressed  in  terms  of  the  time  rate  of 
change  of  the  average  solute  concentration,  q,  of  the  particle.   It  can  be 
shown  as  the  following  equation : 
d  q        kf  Ap 


(  C  -  C5  )  12.3) 


d  t        Vp  p 
where  C  is  the  concentration  of  the  solute  in  the  liauid  expressed  in 
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moles/liter;  Cg  is  the  equilibrium  concentration  in  the  liquid  phase:  kg   is 
the  film  transfer  coefficient;  and  Ap,  Vp,  and  p   are  the  surface  area,  volumn, 
and  density  of  the  particle,  respectively. 

The  average  solute  concentration  of  the  entire  particle  with  a  radius  of 
R(cm>  is  obtained  by  integrating  the  pointwise  concentration  over  the  volume 
of  the  particles,  thereby  yielding 
3 


R3  J 


q  r2  dr  (2.9) 


For  a  spherical  particle,  assuming  symmetry  in  two  directions,  tranport  within 
the  solid  phase  is  given  by  : 

9q        Ds     9         3q 
= (  r2  )  (2.10) 


8t  r2    3r  dr 

or 

8q  a2q      2    3q 

=  D«  (  * )  (2.11) 


at  ar2      r  ar 

q  is  the  dimensionless  solid  phase  concentration  at  an  internal  radius,  r. 
and  Ds  is  the  solid  phase  diffusion  coefficient.  The  mass  balance  for  the 
batch  reactor  is 


d  C  d  q 

V  =   -  W  (2.12) 

d  t  d  t 


V  is  the  volume  of  solution,  and  W  is  toatl  weight  of  carbon.  The  initial  and 

boundary  conditions  for  the  liquid  and  solid  phases  are  given  as 

e       t   =   0    :      C   =   CQ  (2.13) 

@        t    =   0,0<r<R      :    q    =    0  (2.14) 

9        t    >   0    .    r    =   0      :       aq/9r=0  (2.15) 
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t  >  0  ,  r  =  R 


R2  kf  3 

(  C  -  Cs  )  =  |  q  tz   dr  (2.16) 


p  at 

I   r  =  R  ,  Cs  =  f  Cqs)  (2.17) 

Equations  2.13  and  2.14  imply  that  the  initial  concentration  of  solute  is  CG 
in  the  liquid  phase  and  zero  in  the  solid  phase,  respectively.  Equation  2.15 
implies  that  at  any  time  greater  than  zero,  the  concentration  gradient  at  the 
center  of  the  particle  is  zero.  Equation  2. 16  gives  the  mass -balance  on  the 
adsorbent  particles  which  can  be  obtained  from  substituting  Equation  2.9  into 
Equation  2.8.  Equation  2.17  is  the  isotherm  relation  between  the  surface  and 
solution  concentration  of  the  solute.  The  three  parameter  equation.  Equation 
2.6,  is  used  to  describe  equilibrium,  qs,  at  the  external  surface. 
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2-1-3  Fixed  Bed  Model 

In  deriving  the  fixed  bed  model  the  following  assumptions  are  made. 
(1>   Local  equilibrium  exists  at  the  external  surface  of  adsorbent 

particles. 
<2>   Liquid  diffusion  resistance  occurs  at  the  external  surface  and  can  be 

described  as  film  transfer. 

(3)  An  adsorbent  particle  is  a  homogeneous  and  spherically  shaped  solid. 

(4)  The  liquid  phase  is  uniform  across  any  bed  section. 

<5)  The  effect  of  axial  dispersion  is  negligible  in  the  range  of  flowrate 
typical  of  a  liquid  adsorption  system. 

(b)  Diffusion  coefficients  and  rate  parameters  are  independent  of  concen- 
tration. 

Equations  incorporated  in  describing  the  fixed  bed  model  include  Equation 

2.11  and  Equations  2.13  to  2.17  as  solid  phase  conditions.   Other  eauations  in 

the  liquid  phase  conditions  are  given  below: 

dC         dC         1  -  ep 
=  -  v  -  [  3  3  kf  <  C  -  Cm  )  <2.i8) 


at       az       «p  r 

I   t=O.0<Z<LD   :   C=0  (2.19) 

@   t  >  0  ,  Z  =  0       :C=CG  (2.20) 

Equatiion  2. IS  gives  the  mass  balance  of  the  solute  at  any  time  in  the  bed. 
The  rate  of  change  in  the  concentration  of  solute  is  equal  to  the  sum  of  the 
convective  change  in  the  concentration  inside  the  bed  and  rate  of  mass  tran- 
sfer through  the  film  to  the  solid  surface.   Here,  v  is  the  superficial  liquid 
velocity*  «B,  the  porosity  of  bed,  and  Z,    the  bed  depth.   Equation  2.19 
represents  the  initial  condition  to  the  bed  at  tine  zero.   Equation  2.20  reDr- 
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esents  that  at  any  time  equal  to  and  greater  than  zero,  the  solution  concen- 
tration at  the  entrance  to  the  bed  is  equal  to  the  initial  concentration.   Lb 
is  the  length  of  bed,  and  Co  is  initial  solution  concentration. 

The  above  equations  are  transformed  into  non-dimensional  form  by  introdu- 
cing following  variables: 
7  =  r  /  R 
q  =  q  /  qe 
Is  '   <!s  'is 
C  =  C  /  C0 

Cg  ■  Cg  /  Cq 

Z  =  Z  /  Lb 

7  ■  Lb  /  v 

Dg  =  qe  (  1-  eg  )  /C0 

T'  =  t  /  y  Dg 

E<j  ■  Ds  Dg  7  /  r2 

St  =  kfr(l-eg)/R 

The  non-dimensional  fixed  bed  model  then  becomes  : 

3q     Ed     3         9q 

=  — (  r2  )  (2.21) 

9T'    r2   3r         9r 

8  T'  =  0  ,  0  <  r  <  1  :   q  =  0  (2.22) 
@   T'  >  0  ,  r  =  0   :   3q/9r  =   3  (2.23) 

t  T»  >  a  ,  r  »  1   : 

3     P 

St  <  C  -  Ca  >  =  —3-    q  f 2  dr  <2.24) 

8r   ; 

•  r  •  1  :  Cs  =  f(q5) 

9C  8C  _ 

Dg  —  -  3  Dg  St  (  C  -  Cs  )  (2.25) 


9T  az 
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9  T'=0,0<2<1   :  C  =  0  (2.26) 

I   T'>0,Z  =  0:C  =  f(T')  (2.27) 

The  differential  equations  in  the  fixed  bed  model  have  been  sucessfully 
solved  by  the  use  of  orthogonal  collocation  method  (Liu  and  Weber,  1977; 
Crittenden  et  al.,  1980).   The  collocation  method  presents  one  of  the  several 
other  weighted  residual  methods.   In  this  method,  the  unknown  solution  of  a 
differential  equation  is  approximated  by  a  trial  function  having  constants 
and/or  functions.   When  the  trial  function  is  substituted  into  the  differential 
equation,  the  residual  is  forced  to  zero  at  collocation  points.   The  orthogo- 
nal collocation  is  a  special  case  of  the  collocation  method.   The  trial  fun- 
ctions are  a  set  of  orthogonal  polynomials  and  the  collocation  points  are  the 
roots  of  these  polynomials.   By  the  application  of  orthogonal  collocation,  the 
partial  differential  equations  are  reduced  to  first  order  ordinary  differen- 
tial equations. 
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2-2  Multicomponent  Adsorption  Models 

2-2-1  Equilibrium  Models 

The  most  •iHpliUt  model  for  describing  adsorption  equilibrium  in  a  multi- 
component  system  is  the  Langmuir  model  for  competitive  adsorption  (Butler  and 
Ockrent,  1930,  Markam  and  Benton,  1931).   The  basic  assumptions  of  this  model 
is  the  same  as  the  Langmuir  model  for  a  single  solute  system.   The  equation  is 
given  as 

Qi  bi  Cj 
qi  =  (2.28) 

1  ♦  2b]  C] 

The  surface  concentration  of  adsorbate  i,  q±,    is  expressed  as  a  function  of 
the  solution  concentration  of  all  N  adsorbates  in  the  mixture.   The  parameter, 
Qi,  bj.,  are  obtained  directly  from  single  solute  system. 

The  Langmuir  competitive  equation  has  been  modified  by  Schay  et  al.  <19S5) 

by  the  use  of  interaction  terms  to  improve  the  description  of  equilibrium 

data.   The  model  is  represented  as 

Qi  bi  (Ci/^) 

qi  =  (2.29) 

1  +  S  bi  (Ci/i7i) 

where  ^  is  the  interaction  tern  for  each  solute,  this  term  is  evaluated  by 

correcting  the  model  to  multicomponent  experimental  data.   The  main  reason  of 

introducing  interaction  terms  is  that  interaction  between  components  in  a 

mixture  can  reduce  the  adsorption  of  some  components  to  levels  below  that 

predicted  from  the  single  solute  isotherm. 

Mathews  (1975)  has  extended  the  three  parameter  isotherm  to  multicomponent 
adsorption  which  is  expressed  as 
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Ai  Ci 

qi  ■   -■  (2.30) 

1  *  S  Bj  Cj* 

The  values  of  Ai,  Bj  ,and  fii   are  obtained  from  single  solute  system  that  are 
described  by  the  three  parameter  isothera.   Also,  Mathews  has  proposed  the  use 
of  an  interaction  term  ,  as  done  by  Schay  et  al.  (1975),  to  adequately  des- 
cribe a  multicomponent  system.   As  a  result.  Equation  2.30  is  then  given  as 

Ai  (Ci/TTi) 

qi  =  jr  (2.31) 

1  *   Z  Bi  <Cj/i7j)  J 

The  interaction  coefficients,  *7i,  are  calculated  by  the  correlation  of  mixture 

equilibrium  data. 

The  ideal  adsorbed  solution  model  (IAS  model)  has  been  developed  from 
dilute  solution  thermodynamics  by  Radke  and  Prausnitz  (1972a)  using  the  method 
of  Myers  and  Prausnitz  (1965)  for  gas  mixtures.   The  model  predicts  multicom- 
ponent  equilibrium  by  using  single  solute  isotherm.   The  basic  assumptions  of 
the  IAS  models  are: 

(1)  The  available  specific  surface  area  is  identical  for  all  adsorbates. 

(2)  The  adsorbent  is  thermodynamically  inert, 

(3)  Single  solute  isotherm  data  are  required  to  extremely  low  concentra- 
tions, generally  below  5  x  10 ~5M  and  varying  with  the  solutes. 

(4)  The  adsorbed  phase  forms  an  ideal  solution. 

For  strongly  adsorbed  solutes,  predicted  results  are  not  in  good  agreement 
with  experimental  data. 
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2-2-2  Batch  Kinetic  Models 

Weber  and  Morris  (1964),  in  their  studies  on  competitive  adsorption  from 
solution  containing  two  solutes,  found  that  each  solute  adversely  affected  the 
rate  of  adsorption  in  the  presence  of  the  other  solute.   But,  the  overall  rate 
of  adsorption  was  greater  than  those  for  either  of  the  individual  substances 
from  its  pure  solution. 

Tien  and  Thodos  (1976)  described  ion-exchange  kinetics  using  the  homoge- 
neous solid  phase  diffusion  and  a  Freundlich  type  isotherm.   Brecher  et  al. 
(1973)  provided  finite  difference  approximations  for  the  model  using  BET 
isotherm.   Finite  difference  approximations  were  developed  for  any  arbitrary 
isotherm  by  Mathews  and  Weber  (1975).   They  also  extended  the  model  to  multi- 
component  systems  and  verified  the  model  for  batch  adsorption  for  a  wide  range 
of  binary  mixtures.   This  approach  has  been  successfully  applied  to  several 
solutes  and  mixtures  in  batch  reactors  and  fixed  beds  (Crittenden  and  Weber, 
1978;  Lee  et  al,  1980;  Weber  and  Liu,  1980;  Weber  and  Pirbazari,  1981;  Thacker 
et  al,  1981). 
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2-2-3   Fixed  Bed  Models 

Yen  and  Singer  (1984)  applied  the  IAS  model  with  a  modified  calculation  to 
test  ten  sets  of  binary  and  tertiary  phenolic  mixtures.   The  Langmuir  compe- 
titive model  was  used  for  comparison.   The  IAS  model  was  found  to  be  success- 
ful in  precisely  describing  the  competitive  adsorption  equilibria  of  phenolic 
mixture  and  was  proven  to  be  superior  to  the  Langmuir  competitive  model  in 
all  cases  of  studies. 

Crittenden  and  Weber  (1978a,  1978b,  1978c)  have  studied  single  solute  and 
bisolute  adsorption  of  the  same  solutes  and  adsorbent  as  used  by  Mathews 
(1975).   They  used  the  homogeneous  solid  diffusion  model  as  a  basis.   The  film 
transfer  and  the  solid  phase  diffusion  coefficients  were  calculated  from 
single  solute  batch  kinetic  experiments  as  mentioned  in  the  last  section.   The 
model  successfully  described  the  single  solute  fixed  bed  but  not  from  multico- 
mponent  fixed  bed  experiments.   The  discrepancies  observed  were  thought  to  be 
due  to  (1)  experimental  data  scatter;  (2)  poor  equilibrium  description;  (3) 
failure  of  the  assumption  of  independent  diffusion  of  solute. 
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2-3  Factors  Influencing  Adsorption 

2-3-1  Temperature 

The  temperature  at  which  an  adsorption  process  is  conducted  will  affect 
both  the  rate  of  adsorption  and  the  extent  to  which  adsorption  occurs.   Adsor- 
ption rate  parameters  generally  increase  with  increased  temperature  and  dec- 
rease with  decreased  temperature.   The  extent  of  adsorption  will  increase 
with  a  decrease  in  the  lower  temperature  and  decrease  with  an  increase  in  the 
higher  temperature  due  to  being  en  exothermic  nature  of  the  adsorption  reac- 
tion. 

Because  water  is  adsorbed  from  the  surface  when  adsorption  from  aqueous 
solution  occurs,  heat  effects  for  the  process  are  somewhat  smaller  in  compari- 
son with  gas-phase  phase  adsorption.  Heat  of  gas-phase  adsorption  generally  is 
several  Kcal  per  mole.   Thus  small  variations  in  temperatures  do  not  tend  to 
alter  the  adsorption  process  to  a  significant  extent  for  solid  phase  adsor- 
ption. 

Suzuki  and  Kawazoe  (1975)  studied  adsorption  of  fifteen  kinds  of  volatile 
organics  on  activated  carbon  from  aqueous  solution  carried  out  in  a  batch 
system.   They  assumed  that  the  adsorption  rate  was  controlled  by  intraparticle 
diffusion  and  the  solid  phase  diffusion  coefficients  for  each  organic  was 
determined  by  applying  the  concentration-time  curves.  The  solid  phase  diffu- 
sion coefficient,  Ds  (cm^/sec),  thus  obtained,  is  successfully  correlated  to 
the  ratio  of  the  boiling  point  of  adsorbate  and  to  adsorption  temperature  as 

Tb 
Ds  =  a  exp  (-  b  )  (2.32) 
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where  a  and  b  are  constants  from  experimental  data. 

Sudo  and  Suzuki  <1978)  determined  the  concentration  dependence  of  the 
solid  phase  diffusion  coefficient  during  adsorpion  from  the  aqueous  phase. 
The  solid  phase  diffusion  coefficient,  Ds,  was  determined  as  a  function  of  the 
amount  of  solute  adsorbed  on  the  carbon,  q.   Based  on  the  data  obtained,  it 
follows  that  Ds  increased  with  q  and  that  logD6  vs  q  has  a  constant  slope  for 
all  organics  investigated  within  the  range  of  the  amount  adsorbed.   The  follo- 
wing relationship  represented  the  result. 

D5  ■  Dso  exp(a  q)  (2.33) 

where  Dso  was  determined  from  Ds  values  extropolated  to  where  q  =  0. 

Suzuki  and  Fujii  (1982)  studied  concentration  dependence  of  surface  diffu- 
sion coefficients  of  propionic  acid  in  activated  particles.   They  conducted  a 
separate  set  of  experiments  from  which  adsorption  relations  were  determined  by 
batch  adsorption  at  283°K,  293°K,  303°K.  313°K  for  a  wide  ranges  of  concentra- 
tion from  2  x  10"3  to  20  mole/m^.   Strong  dependence  of  Ds  on  q  is  partially 
interpreted  in  terms  of  the  change  of  heat  of  adsorption  with  surface  coverage 
as  determined  from  separate  equilibrium  runs.   The  Ds  values  were  calculated 
by  a  steady  state  technique  for  determining  dependence  of  the  Ds  on  amount  of 
adsorption  which  was  proved  to  be  useful.   The  steady  state  diffusion  experi- 
ment was  performed  with  varying  concentrations  of  propionic  acid  solution 
passed  through  activated  carbon  pellets.   The  following  equation  represents 
the  relation  between  Ds,  T  and  q. 

Qo 

Ds  =  Dso  a*P  1  "  1"  <a  q>3  (2.34) 

R  T 

Qo  is  obtained  from  the  following  equation. 

Qst,a  =  -  Qo  ln(a  q)  (2.35) 
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where  Ost.a  is  isosteric  heat  of  adsorption.   It  can  be  defined  as 

d  In  C 

Ost.a  ■  8  <2-36> 

d  (1/T) 
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2-3-2  Other  factors 
Surface  Area 

Adsorption  is  a  surface  phenomenon  as  such,  the  extent  of  adsorption  is 
proportional  to  specific  surface  area.   The  definition  of  specific  surface  area 
is  the  portion  of  surface  area  that  is  available  for  adsorption.   Thus  the 
amount  of  adsorption  accomplished  per  unit  weight  of  a  solid  adsorbent  is 
greater  the  more  finely  divided  and  the  more  porous  solid. 

Agitation 

The  rate  of  adsorption  is  controlled  by  either  film  transfer  or  surface 
diffusion,  depending  upon  the  intensity  of  agitation  in  the  system.   If  rela- 
tively little  agitation  occurs,  the  film  transfer  will  likely  be  the  rate- 
limiting  step  because  of  the  surface  film  around  the  particle  will  be  thick. 
For  a  continuous-flow  system  of  flowrate  less  than  10  gal/ft~-min,  the  film 
transfer  will  be  rate-limiting.   If  adequate  mixing  is  provided  like  in  batch- 
type  contacting  system,  surface  diffusion  is  generally  rate  limiting. 

Size  of  Carbon  Particles 

Morris  and  Weber  (1964,1964a)  reported  that  the  rate  of  adsorption  corre- 
lated with  the  inverse  of  the  square  of  the  diameter  of  carbon  particles. 
Specifically,  adsorption  rates  increase  with  as  particle  sizes  decrease.   Thev 
also  observed  a  significant  increase  in  adsorption  capacity  with  decrease  in 
carbon  partical  size. 

In  the  column  system,  the  time  required  to  reach  the  breakpoint  increased 
linearly  with  decrease  in  particle  size.   After  the  breakpoint,  however,  the 
approach  to  the  exhaustion  capacity  of  the  carbon  column  speeded  up  with 
decreasing  particle  size  (Martin  1978). 
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Solubility  of  Adsorbatea 

In  general,  an  inverse  relationship  can  be  anticipated  between  the  extent 
of  adsorption  of  a  solute  and  its  solubility  in  the  solvent  from  which  adsor- 
ption occurs.   The  greater  the  solubility,  the  stronger  the  solute-solvent 
bond  and  the  smaller  the  extent  of  adsorption. 

EH 

Adsorption  of  typical  organic  pollutants  from  water  is  increased  with 
decreasing  pH.   In  many  cases  this  may  result  from  neutralization  of  negative 
charge  at  the  surface  of  the  carbon  with  increasing  hydrogen-ion  concentra- 
tion, thereby,  reducing  hindrance  to  diffusion  and  making  available  more  of 
the  activated  surface  of  the  carbon. 


CHAPTER  3. 
EXPERIMENTAL  MATERIALS  AND  METHODS 

3-1  Materials 

3-1-1  Adsorbents 

The  adsorbent  used  in  this  research  was  granular  activated  carbon.   The 
physical  properties  of  the  carbon  are  listed  in  Table  3.1. 

Filtrasorb-400  were  chosen  for  all  experiments  in  this  study.   The  carbon 
is  received  as  12/40  (particles  passing  through  a  U.S.  Standard  Mesh  Size 
No. 12  screen  but  retained  on  No. 40  screen)  and  is  seived  to  20/25  to  give  a 
narrower  size  range  for  batch  kinetic  and  fixed  bed  studies.   The  geometric 
mean  diameter  of  a  particle  is  given  by  <d]>d2>ly'2  where  d\    is  the  mesh  size  of 
upper  sieve  and  d-?   is  the  mesh  size  of  the  lower  sieve.   Based  on  this, 
particles  passing  through  U.S.  mesh  size  No. 20  and  retained  on  U.S.  mesh  size 
No. 25  had  a  diameter  of  772  microns.   For  equilibrium  studies,  some  carbon 
particles  were  ground  to  provide  a  small  particle  size  200/230.   Before  expe- 
riments, the  carbon  was  washed  several  times  in  distilled  deionized  water, 
dried  to  constant  weight  in  an  oven  at  a  temperature  of  103°C  and  kept  at  room 
temperature  in  airtight  containers  pending  use. 
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Table  3.1 


Properties  of  Activated  Carbon 


Manufacturer 

Calgon  Co.  PA 

U.S.  Mesh  Size 

20  /  25 

Raw  Material 

Bituminous  Coal 

Physical  Properties 

surface  area,  »2/ga 

1000  -  1100 

apparent  density,  gm/cc 

2.1 

particle  density,  gm/cc 

1.35 

wetted  in  water 

effective  size,  mx 

0.71  -  0.84 

uniformity  coefficient 

1.9  or  less 

bachwashed  and  drained 
density,    lb/ft3 
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3-1-2  Adsorbates  :  Phenol  and  Parachlorophsnol  (PCP) 

The  adsorbates  used  in  this  research  were  phenol  and  parachlorophenol 
(PCP)  with  their  properties  listed  in  Table  3.2.   Both  adsorbates  were  sup- 
plied by  J.T.  Baker  Chemical  Co.  in  crystal  for»  stored  in  dark  glass  bottles. 

All  stock  solutions  were  prepared  with  organic-free  distilled  water.   The 
system  to  obtain  this  water  was  a  4-stage  filter  which  is  able  to  remove 
organic  compounds  in  tap  water  and  produce  ultra  pure  distilled  water.   The 
stock  solutions  were  stored  in  a  freezer  and  diluted  to  the  desired  concentra- 
tion prior  to  use.   The  pH  values  through  experiaents  were  kept  constant  at  7. 


Table  3.2 


Adsorbates  and  Their  Properties 
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Phenol 

PCP 

Supplier 

J.T.  Baker 

J.T.  Baker 

Commercial  Co. 

Commercial  Co. 

Grade 

Baker  Reagent 

Grade 

Baker  TO 

Foraula 

C6H50H 

C6H4OHCI 

M.W. 

94.1 

128.58 

pKa 

9.90 

9.18 

Melting  point 

43°C 

42-44°C 

Boiling  point 

181.75°C 

- 

• 

Diffusivity 

9.12  at  20°C 

8.71  at  20°C 

in  water 

<10-6  cm2/sec) 

10.43  at  25°C 

9.96  at  25°C 

11.83  at  30°C 

11.32  at  30°C 

13.33  at  35°C 

12.74  at  35°C 
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3-2  Methods 

3-2-1   Analytical  Methods 

The  concentrations  of  phenol  and  PCP  solutions  were  analyzed  by  ultravio- 
let absorption  spectrophotometry  using  a  Baush  &  Lomb  Spectronic  710  Spectro- 
tometer.   The  wavelengths  of  «oxi«u>  absorbance  for  phenol  and  PCP  were  deter- 
mined in  preliminary  tests.   Several  standard  solutions  of  phenol  and  PCP  were 
prepared  after  the  wavelengths  of  maximum  absorbance  were  Measured.   To  obtain 
significant  absorbance  when  analyzing  low  concentration  samples,  50  mm  photo- 
cells were  used,  whereas  more  concentrated  samples  were  analyzed  with  a  10  mm 
photocell.   Calibration  curves  of  standard  solutions  were  determined  for  light 
path  lengths  of  10  mm  and  SO  ma. 

Phenol  was  analyzed  at  a  wavelength  of  268  nm  and  PCP  at  279  nm .   Table 
3.3  shows  the  selected  wavelengths  and  the  corresponding  molar  absorptives 
which  are  determined  from  the  slopes  of  the  calibration  curves. 

Mixtures  of  phenol  and  PCP  were  analyzed  at  wavelengths  of  268  nm  and  279 
nm.  The  procedure  developed  by  Friedel  (1951)  was  employed  to  compute  indivi- 
dual concentrations  of  phenol  and  PCP.  Note  that  underlying  assumption  of  the 
Friedel 's  procedure  is  that  the  absorbance  of  a  mixture  can  be  expressed  as  a 
linear  combination  of  that  of  it  compounds  at  a  given  wavelength.  As  a 
result,  the  following  set  of  simultaneous  equations  are  used  to  determine  the 
concentration  of  phenol  and  PCP  in  a  mixed  system. 

*268  =  «P1  Cp  *  epCPl  Cpcp  (3.1) 

A279  =  ep2  Cp  ♦  epcP2  Cpcp  (3.2) 

where 
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A268'A279  =  the  absorbance  of  the  mixed  solution  at  wavelengths  of 

268,  279  nm,  respectively. 
Cp,  Cpcp  ■  concentration  of  phenol  and  PCP.  respectively. 
<P1**PCP1  =  Molar  absorptives  of  phenol  and  PCP  at  a  wavelength 

of  268  nm. 
GP2»SPCP2  =  Molar  absorptivities  of  phenol  and  PCP  at  a  wavelength 
of  279  nm. 
The  concentrations  of  Phenol  and  PCP,  i.e.,  Cp  and  Cpcp.  respectively,  can  be 
determined  by  solving  the  above  simultaneous  equations. 
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Table  3.3 


Molar  Abaorptivities  of  Adsorbates  (1/nole-cm) 


wavelength  (n»)  phenol  PCP 


268  1463  887 


279  592  1438 
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3-2-2  Equilibrium  Experiments 

Equilibrium  experiments  were  carried  out  by  the  standard  bottle-point 
method.   Different  amounts  of  particles  (200/230)  were  weighed  carefully  and 
put  into  a  series  of  125  ml  glass  bottles  which  contained  100  ml  of  phenol 
and/or  PCP  solution  with  concentration.  Carbon  was  prepared  by  washing  samples 
several  times  with  deionized  water  to  remove  fines.   Glass  bottles  were  washed 
throughly  and  dried  in  an  oven  at  103°C.  Carbon  dosages  varying  from  0.1  -  2.2 
gm  were  added  to  the  bottles  prepared  earlier  as  a  adsorbate.   The  bottles 
were  sealed  with  caps  to  ensure  airtight  and  placed  in  a  Controlled  Environ- 
ment Incubator  Shaker  (Labline  Instruments  Inc.,  Ill)  to  provide  agitation  at 
desired  temperatures. 

In  each  experimental  run,  the  initial  solute  concentration,  Co,  was  recor- 
ded at  the  outset  of  the  run,  followed  by  measurement  of  the  bulk  solute 
concentration,  Ce,  on  a  daily  basis;  the  amount  of  solute  adsorbed  per  unit 
mass  of  activated  carbon,  q,  was  thus  determined  as 

Cq  -  Cg 

qe  s   V  (3.3) 

W 

where  W  is  the  total  weight  of  the  carbon,  and  V  is  the  volume  of  solution. 
It  is  worthwhile  mentioning  here  that  measurement  of  solute  concentration 
involved  a  two-step  process.   In  step  one,  the  samples  were  filtered  through  a 
23  ftM   membrane  filter  (Gelman  Science  Co.)  to  remove  any  suspended  carbon;  the 
solute  (filtrate)  concentrtation  was  then  measured  subsequently.   The  amount 
of  each  solute  adsorbed  on  the  carbon  at  equilibrium  was  calculated  from  the 
mass-balance  Equation  3.3.   Each  bottle  represents  one  point  on  the  equilib- 
rium curve.   By  using  different  initial  concentrations  and  by  varying  the 
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carbon  dosage,  the  isotherm  profile  can  be  determined.   The  isotherm  data  are 
plotted  typically  on  log-log  paper.   The  concentration  remaining  in  solution 
after  equilibrium,  Co,  are  plotted  on  the  abscissa  ;  the  amount  of  solute 
adsorbed  per  unit  weight  of  activated  carbon,  qe,  are  plotted  on  the  ordinate. 
The  resulting  plots  represent  the  adsorption  isotherm.   A  series  of  adsorption 
isotherm  data  were  then  determined. 

To  provide  an  environment  with  constant  temperature  for  equilibrium 
experiments,  the  Controlled  Environment  Incubator  Shaker  was  set  to  a  desired 
temperature.   The  shaker  was  adjusted  so  that  the  adsorbent  was  always  kept 
suspended,  care  being  taken  to  see  that  the  adsorbent  did  not  stick  to  the 
sides  of  the  bottles.   Before  filtration,  the  bottles  were  removed  from  the 
shaker,  and  they  were  partially  immersed  in  a  bath  in  which  the  temperature 
was  maintained  at  the  same  temperature  as  that  of  the  Environment  Incubator 
Shaker. 

Four  different  temperatures,  20°,  25°,  30°.  and  35°C,  were  tested  for 
experiments  with  single  solute.   In  cases  of  mixtures,  the  experiments  were 
conducted  under  only  two  temperatures,  namely,  20°  and  35°C. 
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3-2-3  Batch  Kinetic  Experiments 

Carbon  adsorption  processes  investigated  in  this  study  were  operated  on 
both  batch  and  fixed  bed  basis.   The  former  is  discussed  in  this  section  while 
the  latter  will  be  elaborated  in  next  section. 

Batch  experiments  were  conducted  to  evaluate  the  film  transfer  coefficient 
<kf>  and  the  solid  phase  diffusion  coefficient  <DS)  on  carbon  of  phenol  and 
PCP  at  various  temperatures.   Batch  kinetic  experiments  were  conducted  in  a 
rectangular  vessel  made  of  plexiglass.   The  dimensions  of  the  vessel  were  34 
cm  x  30  cm  x  36  cm  which  are  L,  B,  and  H,  respectively.   Wall  thickness  was 
0.68  cm.   A  steel  impeller  was  located  at  the  center  of  the  vessel  and  7  cm 
from  the  bottom.   A  1000  watt  heater  and  a  temperature-controller  were  used  to 
regulate  the  temperature  of  the  solution.   The  schematic  of  experimental 
facilities  are  illustrated  in  Figure  3.1.   The  rate  experiments  were  normally 
conducted  for  a  duration  of  3  or  3.5  hours. 

A  solution  of  24  liters  was  prepared  from  tap  water  passed  through  a  three 
carbon  columns  to  remove  suspended  and  organic  matter.   The  depth  of  each 
solution  from  the  bottom  was  about  23.5  cm.   In  no  experiments  were  crystals 
of  solute  directly  introduced  into  the  vessel.   In  all  cases,  stock  solutions 
with  the  desired  concentration  were  made  up.   The  impeller  was  maintained  at  a 
speed  of  700  rpm  to  agitate  the  solution  to  reach  equilibrium  within  the 
vessel.   A  fixed  amount  of  solution  was  withdrawn  to  determine  the  initial 
concentration  after  the  impeller  ran  steadily.   A  predetermined  amount  of 
carbon  was  poured  abruptly  into  the  reaction  vessel  at  which  moment  the  time 
started  counting.   Samples  were  withdrawn  from  the  vessel  with  a  pipet  to 
analyze  concentration  at  selected  time  intervals. 
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The  size  of  carbon  particle  used  in  experiments  ranged  fro*  U.S.  mesh  size 
No. 20  to  No. 25.   The  total  weight  of  the  carbon  used  was  6  gm.   The  solution 
was  preheated,  if  necessary,  to  the  desired  temperature.   A  temperature  con- 
troller was  used  to  keep  the  temperature  at  a  constant  value.   For  a  step 
change  in  temperature,  solution  was  held  constant  in  the  first  hour  followed 
by  the  use  of  a  heater  or  dry  ice  for  an  increase  or  decrease  in  temperature. 
It  took  15  minutes  to  raise  the  temperature  from  20°C  to  35°C.   To  decrease 
temperature,  approximate  1.5  lb  of  dry  ice  was  introduced  into  the  vessel,  and 
it  took  20  minutes  to  bring  down  the  temperature  from  35°  to  20°C.  For  con- 
stant temperatures,  four  runs,  each  with  temperature  of  20°,  25°,  30°,  and 
35°C,  respectively,  were  carried  out  for  single  solute  studies;  two  runs,  each 
with  temperature  of  20°  and  35°C,  respectively,  were  for  mixture.   For  step 
changes,  two  runs  were  done.   One  with  the  temperature  increased  from  20°  to 
35°C  and  the  other  with  temperature  decreased  from  35°  to  20°C. 
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3-2-4  Fixed  Bed  Experiments 

Figure  3.2  depicts  the  schematic  of  the  fixed  bed  experiment.   The  fixed 
bed  was  a  plexiglass  tube  with  an  inner  diameter  of  3.1  cm  and  a  height  of  6 
feet.   Glass  beads  with  2  mm  diameter  were  filled  at  the  bottom  of  the  column 
as  a  velocity  distributor  and  maintained  uniform  upflow  velocity.    Its  fun- 
ction was  to  distribute  the  solution  if  the  fixed  bed  employed  upflow  of 
liquid.   In  the  research,  all  experiments  were  conducted  with  downflow.   A 
140-mesh  brass  screen  separated  the  carbon  in  the  column  from  the  glass  beads 
in  the  velocity  distributor.   A  30-gallon  tank  supplied  a  continuous  flow  of 
water  to  the  top  of  the  column  as  downflow  of  solution;  the  water  was  produced 
by  having  tap  water  run  through  a  three-foot  long  filter  filled  with  activated 
carbon  to  remove  organics  and  other  adsorbates  from  water.   Two  varible  flow 
pumps  were  manipulated  to  force  the  water  as  well  as  the  concentrated  solution 
at  a  combined  constant  flowrate  of  125  ml/»in.   The  solution  stored  in  a 
solution  jar  was  prepared  at  a  higher  concentration  than  -that  was  required. 
Simultaneous  adjustments  on  the  flows  of  the  concentrated  solution,  and  dis- 
tilled water  were  performed  to  obtain  a  desired  influent  solution  concentra- 
tion.  Constant  checks  on  this  influent  concentration  were  required  to  ensure 
a  minimum  variation  in  flowrate.   The  fixed  bed  experiments  were  generally 
conducted  until  the  ratio  of  effluent  concentration  to  influent  concentration 
reached  0.95.   The  durations  of  the  experiment  were  usually  from  10  to  15 
hours. 

The  desired  influent  solution  temperature  was  obtained  by  mixing  the 
concentrated  solution  at  an  ambient  temperature  with  the  pre-adsorbed  tap 
water  at  a  higher  temperature;   the  determination  of  the  distilled  water 
temperature  involved  trial-and-error  efforts.   As  an  example,  distilled  water 
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was  preheated  to  a  temperature  of  54°C  so  as  to  produce  a  resultant  solution 
with  a  desired  temperature  of  35°C.   Prior  to  the  experiment  the  bed  was 
fluidized  for  30  minutes  to  remove  trapped  bubbles  in  the  actived  carbon  zone. 

For  constant  temperature  studies,  temperatures  were  selected  at  20°  and 
35°C  for  runs  with  single  solution.   For  step  change  in  temperature  studies, 
the  temperature  was  held  constant  for  the  first  three  hours  of  the  experiment 
before  the  step  change.   It  took  30  minutes  to  raise  the  temperature  from  20° 
to  35°  and  25  minutes  to  reduce  the  temperature  from  35°C  to  20°C.  Temperature 
settings  were  the  same  as  those  described  in  the  kinetic  experiment  section. 

Effluent  solutions  were  collected  and  analyzed  for  determination  of  concen- 
tration.  Breakthrough  curves  were  hence  obtained  by  plotting  the  ratio  of 
effluent  to  influent  concentration,  Ci/Co,  versus  time  ,T. 
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CHAPTER   4. 
SINGLE  SOLUTE  ADSORPTION 

4-1   Eauilibrium  Studies 

The  three  oarameter  equation,  which  reduces  to  the  linear,  or   Langnuir 
isotherm ,  was  employed  to  correlate  the  equilibrium  relationship  oetween  the 
surface  concentration,  q«,  and  the  bulk  solution  concentration,  Ce.   The 
equation  mentioned  in  Section  2-1  is 

A  x  Ce 

Qe  e  (2. 6 J 

1  *  B  x  Ce£ 

The  parameters.  A,  S,  and  &,    which  were  determined  by  the  best  statistical 

representation  of  the  experimental  data,  could  be  estimated  by  a  nonlinear 

parameter  estimation  technique  developed  by  Mathews  (1975). 

Experimental  results  of  qe  and  Ce  obtained  in  Section  3-2-2  were  fitted  bv 
Equation  2.6  with  parameter  values  summarized  in  Table  4.1  at  the  four  temcera- 
tures,  20°,  25°,  30°,  and  35°C.   The  corresponding  data  are  plotted  in  Figures 
4.1  through  4.8.   In  these  figures,  it  is  observed  that  the  extent  of  adsor- 
ption increases  with  a  decrease  in  the  temperature:  this  phenomenon  conforms 
the  well-known  fact  that  adsorption  is  an  exothermic  process  in  general. 

Since  the  three  parameter  equation  is  suitable  only  for  a  fixed  temoera- 
ture,  its  applications  are  of  limited  use  in  practice.   As  a  remedy,  an 
attempt  was  made  to  find  the  temperature  dependencies  of  A.  3,  and  &   ,m 
Equation  4.1.   The  method  is  similar  to  that  proDosed  by  Danner  (1933),   The 
Five  parameters.  0o,  S0,  &>.    Kl,  and   K2  were  evaluated  by  running  a  orosram 
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with  entire  equilibrium  experimental  data.   By  many  trial  steps.  Equation  4.1 

was  found  having  the  closest  correlation. 

0o  exp<Kl/T)  Ce 

qe  ■  —  (4.1) 

1  *  B0  exp<K2/T)  Cep' 

It  s  apparent  that  A,  B,  and  fi   can  be  represented  by  Equation  4.2  through  4.4 
by  applying  T  and  the  five  parameters  as  variables. 

A  =  0o  exp(Kl/T)  (4.2) 

B  =  B0  exp(K2/T)  (4.3) 

fi   ■  /So  (4.4) 

The  resultant  parameter  values  of  0o,  B0.  iSo.  Kl,  and  K2  are  listed  in  Table 
4.2.  Table  4.3  contains  the  values  of  A,  B,  and  fi  calculated  from  Equations 
4.2  to  4.4.  Note  that  in  the  parameter  estimation,  fi  is  maintained  tempera- 
ture independent. 

In  Section  2-3-1,  the  isosteric  heat  of  adsorption,  Qst.a  i  can  De 

calculated  from  Equation  2.36. 

d  In  C 

Qst.a  ■  "  (2.36) 

d  (1/T) 

The  plots  of  InC  versus  1/T  for  q  «  0.2,  0.4,  0.6,  0.8,  and  1.0  mmole/gm  are 

showed  in  Figure  4.3.   From  the  slopes,  the  values  of  Qst'a  were  determined. 

Figure  4.10  illustrated  0st.a  as  a  function  of  q.   From  the  experimental 

results,  the  relationship  can  be  represented  by  Equations  4.5  and  4.6. 

phenol : 

Ost.a  ■   -10625  x  ln(0.167  q)  (4.5) 

PCP: 

Qst.a  ■   -15845  x  ln(0.331  q)  (4.6) 


Table  4.1 


Isotherm  Parameters  for  Phenol  and  PCP 


adsorbate    temp(°C> 


phenol 


20 
25 
30 
35 


36.37 

20.34 

0.7705 

32.53 

19.26 

0.7995 

27.09 

17.64 

0.8069 

22.50 

15.43 

0.8272 

42.23 

24.72 

0.8791 

39.56 

26.89 

0.8654 

15.78 

28.65 

0.8747 

33.16 

29.32 

0.8503 

PCP 


20 
25 
30 
35 


unit  of  Equilibrium  cone,  as  mmole/1 
unit  of  surface  cone,  as  mmole/qm 


Table  4.2 

Adsorption  Equilibrium  Constants  for 
a  Range  of  Temperature 
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Constants 


phenol 


PCP 


Oo 
Bo 

Kl 
K2 


0.1515 

1 . 1090 

15.20 

25.12 

0.7838 

0.8500 

1595. 

1057. 

63.32 

22.31 

Table  4.3 

Isotherm  Parameters  for  Phenol  and  PCP 

Calculated  from  Equations  4.2  through  4.4 


adsorbate    temp(°C) 
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phenol 


20 
25 
30 
35 


35.04 

18.87 

0.7838 

31.98 

18.80 

0.7838 

29.28 

18.73 

0.7838 

26.88 

18.67 

0.7838 

40.89 

27.11 

0.8500 

38.49 

27.07 

0.8500 

36.30 

27.04 

0.8500 

34.30 

27.00 

0.8500 

PCP 


20 
25 
30 
35 


unit  of  Equilibrium  cone,  as  mmole/1 
unit  of  Surface  cone,  as  mmole/g» 
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4-2  Batch  Kinetic  Studies 

The  objectives  of  batch  kinetic  studies  were  to  determine  the  solid  phase 
diffusion  coefficients,  D5,  and  the  film  transfer  coefficients,  kf,  of  phenol 
and  PCP.   Carbon  size  of  20/25,  stirring  speed  of  700  rpm,  reaction  volume  of 
24  liters  and  total  carbon  weight  of  6  qm  were  used  in  all  batch  kinetic 
experiments.   The  initial  concentration  of  phenol  and  PCP  was  2.5  x  10_<i  H. 
The  pH  values  were  kept  at  7.   The  reason  for  keeping  the  pH  values  at  7  is 
that  the  solution  will  contain  almost  iOOfc  of  phneol .   The  pKa  value  of  phenol 
is  9.9.   If  the  pH  value  greatly  less  than  pKa,  the  concentration  of  CijHsO" 
approach  zero. 

The  homogeneous  solid  phase  diffusion  model  proposed  by  Rosen  (1952)  was 

provided  to  predict  the  adsorption  rates  in  batch  reactors.   The  principal 

parameters  required  for  use  of  the  model  are  the  film  transfer  coefficient. 

kf,  the  solid  phase  diffusion  coefficient,  Ds,  and  the  isotherm  parameters. 

The  film  trasfer  coefficients,  kf,  were  determined  by  using  the  phenol  or  PCP 

concentration  data  obtained  during  the  initial  concentration.   Equation  4.7 

shows  the  calculation  of  kf. 

-  ln(C/Co)    (1  -  e)  p         R 

kf  = (4.7) 

t       3  W  /  V      2 

The  values  of  ln(C/Co),  C  as  adsorbate  concentration  at  time  t,  and  C0  as 
initial  concentration,  were  plotted  versus  time  t.   The  linear  portion  of  time 
points  were  fitted  with  a  straight  line.   The  values  of  kf.  thus,  were  calcu- 
lated from  the  slope.   At  longer  contact  time,  the  line  become  nonlinear 
because  of  increasing  influence  of  intraparticie  resistance. 

A  parameter  estimation  procedure  used  in  the  homogeneous  solid  phase  model 
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was  demonstrated  to  provide  accurate  estimates  of  the  values  of  kf  and  D3. 
Intraparticle  diffusion  coefficients  of  the  components  were  estimated  by  mat- 
ching the  predictions  of  the  model  to  the  results  of  batch  kinetic  experi- 
ments.  A  computer  program  written  by  Mathews  (1975)  was  used  to  estimate  Kf 
and  Ds  from  batch  reactor  data.   The  input  data  of  the  program  include  isot- 
herm parameters,  the  volume  of  solution,  total  weight  of  carbon,  diameter  of 
particle,  initial  concentration  of  solution,  and  initial  estimates  for  D5  and 
kf.   This  program  uses  the  principal  axis  (Brent  1971)  for  finding  the  minimum 
of  a  function,  and  search  for  parameter  values  or  kf  and  D«  that  will  minimize 
the  sum  of  the  square  of  the  difference  between  experimental  and  computed 
values  of  adsorption  rates. 

Figures  4.11  through  4.14  show  the  results  for  constant  temperature  expe- 
riments.  In  these  figures  the  solid  lines  represent  the  best  fit  of  the  batch 
reactor  model  to  the  experimental  data.   The  values  of  the  film  transfer 
coefficient,  kf,  and  the  solid  phase  diffusion  coefficient,  D5.  are  listed  in 
Table  4,4.   Comparison  of  the  curves  shows  that  PCP  has  higher  adsorption 
rates  than  those  of  phenol.   From  experimental  results,  both  the  film  transfer 
and  solid  phase  diffusion  coefficients  increase  with  increasing  temperature. 

For  step  change  experiments,  two  tests  were  run  for  phenol  and  PCP,  that 
is,  from  20°  to  35°C  and  from  35°  to  20°C.   The  computer  program  origina- 
lly could  only  predict  the  coefficients  at  constant  temperatures.   A  modifica- 
tion developed  by  Mathews  (1984)  takes  into  account  temperature  as  a  variable 
in  fitting  step  change  experiments.   In  the  modification.  Equation  2.6  was 
first  subsituted  by  Equation  2.1  to  include  temperature  effect  on  the  adsor- 
ption capacity.   The  film  transfer  coefficient,  kf,  and  the  soiid  phase  diffu- 
sion coefficient,  Ds  ,  following  the  Arrhenius'  Law,  were  represented  as 
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Equation  4.8  through  4.11. 

phenol:    kf  =  3.58  x  Q-1975.3/T  <4.8> 

Ds  =  6.4637  x  1014  x  q-15437/T  (4#9) 

PCP:       kf  =  71.144  H  Q-2863.3/T  C4.10) 

D5  =  1.857  x  106  x  e"9473.3/T  (4.11) 

where  T  represents  the  absolute  temperature  °K.   The  above  equations  result 
from  least-square  fitting  of  the  data  in  Table  4.6.   The  values  calculated 
from  the  equations  above  are  listed  in  Table  4.7.   The  plots  of  -in<kf)  and 
-ln(Ds>  versus  1/T  are  shown  in  Figure  4.15. 

It  was  mentioned  that  15  minutes  was  needed  to  raise  the  temperature  from 
20°C  to  35°C  for  the  batch  kinetic  experiments.   In  predicting  the  adsorption 
rates  for  step  changes,  the  program  must  consider  the  time  required  for  the 
step  changes  in  temperature.   After  the  first  hour  of  constant  temperature 
operation,  the  temperature  was  assumed  to  be  25°C  from  60  to  68  minutes.  30°C 
from  68  to  75  minutes  and  35°C  thereafter.   Consequently,  during  each  tempe- 
rature range,  the  isotherm  parameters,  kf  and  Ds,  corresponding  to  that  temoe- 
rature  were  used. 

Figures  4.16  and  4.17  show  the  experimental  data  and  predicted  curves  for 
the  step  change  in  temperature  of  phenol  and  PCP.   For  the  runs  from  20°C  to 
35°C,  the  experimental  data  are  in  good  agreement  with  the  predicted  curves. 
For  the  runs  from  35°C  to  20°C,  the  experimental  data  gives  a  fair  fit  with 
the  predicted  curves.   This  could  be  due  to  the  changes  of  adsorption  rates 
and  stirring  speed  during  the  sublimation  when  putting  dry  ice  in  the  solution 
to  decrease  the  temperature.   The  two  curves  in  Figure  4,14  and  4.15  have 
intersections.   The  intersections  are  based  on  that  adsorption  rates  chanqed 
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abroutly  during  the  temperature  of  solution  dropped  to  lower  temperature  or 
raised  to  higher  temperature. 
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Table  4.4 

Solid  phase  Diffusion  Coefficients  and  Film 

Transfer  Coefficients  for  Phenol  and  PCP 

adsorbate      tamp<°C)        Ds  (c»2/sec)  kf  (cn/sec) 

20  1.100  x  10"8  4.273  x  10"3 

25  2.888  x  10"8  4.634  x  10"3 

phenol 

30  7.000  x  10-8  5.331  x  10"3 

35  1.300  x  10"7  5.818  x  10-3 


20  1.521  x  10"8  4.602  x  10"3 

25  2.214  x  10-8  5.410  x  10'3 

30  4.478  x  10"8  6.326  x  10"3 

35  7.022  x  10"8  7.361  x  10"3 


PCP 
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Table  4.5 

Solid  phase  Diffusion  Coefficients  and  Film  Transfer 

Coefficients  Calculated  from  Equation  4.8  to  4.11 

adsorbate    te»p(°C)          Qs  (cm2/sec>  kf  <cm/sec) 

20              1.155  x  10"8  4.227  x  10"3 

25               2.782  x   10"8  4.733  x  10"3 

30               6.507  x  10"8  5.280  x  10"3 

35               1.481  x  10"7  5.869  x  10"3 

20               1.687  x  10"8  4.602  x  10"3 

25               2.902  x  10"8  5.409  x  10"3 

30               4.904  x  10"8  6.326  x  10"3 

35               8.147  x  10"8  7.362  x  10"3 


phenol 


PCP 


63 


Jo/o 


:,., 


O  r-< 

ft    01 

~ 

4-i      0 

«    E 

^    1 

0 

0 

<u    ' 

CO 

Oh 

rr 

Ma 

o 

o  m 

T 

' 

2 

4J 

ffl 

H 

C    ^ 

a  o 

u  m 

:q 


'0/3 


67 


<b  © 


SE 


Jo/o 


68 


o 

a. 

r\i 

CO 

u    - 

tr 

5 

b  o 

n 

o  m 

X 

on  ^ 

UJ 

u 

2 

CO 

Ss 


J3/0 


69 


(  *M>"I- 


70 


'"  u 

o 

-- 
-I 

rC 

1 1 

0J 

u 

ZZ 

c 

nj  » 

■z 

JZ 

uu 

w 

o i    o 

a. 

rsJrM 

-  «"-, 

"n 

o   = 

O 

a;   u© 

to 

—  ^ 

K 

_J 

s-   c 

o 

0    ra  u 

I 

'SiO     rn 
0J   u"l 

LU 

±J   n  Q 

2 

!0         u 

OS     O 

H 

~          0 

o  cj0 

C  CM 

—  *c 


J0/0 


71 


1 


&  «  o 


n 

u>^ 

CM 

0          -1 

rr 

"-  s 

3 
O 

a.   £© 

T 

o. 

1 — ■ 

S-j     C 

UJ 
2 

0    a  u 

Wo      f> 

1- 

OJ    lO 

—  wo 


kn/o 


72 


4-3  Fixed  Bed  Studies  :  Column 

The  experimental  data  obtained  in  the  fixed  bed  tests  were  applied  to 
judge  the  efficiency  of  the  homogeneous  solid  phase  diffusion  model  in  predic- 
ting breakthrough  curves.   The  operating  data  for  the  fixed  bed  experiments 
are  listed  from  Tables  4.6  through  4.9.   The  input  data  for  the  fixed-bed 
program  for  predicting  breakthrough  curves  include  the  isotherm  parameters,  kf 
and  Ds#  total  weight  of  carbon,  diameter  of  particles,  diameter  of  column, 
total  flow  rate,  surface  loading,  influent  concentration,  and  total  bed 
height. 

Breakthrough  curves  showing  the  plots  of  C/Co  versus  time  indicate  the 
increase  in  the  ratio  of  effluent  to  influent  as  adsorption  in  the  column 
proceeds.   The  general  pattern  of  the  adsorption  breakthrough  curves,  which 
exhibit  like  a  characteristic  S  shape,  are  as  expected  for  all  fixed  bed  runs. 
After  the  trace  of  solute  effluent  was  detected,  the  effluent  concentration 
rose  first  gradually  and  then  sharply  till  it  approached  the  values  of  the 
influent  concentration. 

Figures  4.18  through  4.21  present  experimental  breakthrough  data  and  the 
corresponding  model  prediction  for  phenol  and  PCP  at  temperatures  of  20°C  and 
35°C.   Comparison  of  the  experimental  and  model  predictions  for  these  cases 
indicate  that  the  model  provides  reasonable  projections  of  the  performance  of 
fixed  bed.   From  the  breakthrough  standpoint,  it  is  better  to  operate  the 
column  at  20°C.   Comparing  C/Co  at  0.05,  we  see  that  the  time  for  reaching 
C/Co  at  0.05  is  235  minutes  at  20°C  and  220  minutes  at  35°C  for  phenol  as  weii 
as  for  PCP  460  minutes  at  20°C  and  430  minutes  at  35°C.   The  time  for  C/Co 
reaching  0.10  is  270  minutes  at  20°C  and  250  minutes  at  35°C  for  Dhenol  and 
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510  minutes  at  20°C  and  450  minutes  at  35°C  for  PCP. 

The  Adsorption  capacity  was  higher  at  lower  temperature  from  equilibrium 
experiments.   But  the  adsorption  rates  increased  with  increased  temoerature  as 
seen  from  batch  kinetic  experiments.  These  factors  affecting  adsorption  give 
adverse  influence  on  breakthrough  curves.   From  the  figures,  it  can  be  seen 
that  C/Co  values  at  a  temperature  of  20°C  are  lower  than  those  at  35°C  at 
first  half  profile;  however,  the  values  are  higher  after  that  moment. 

Step  change  experiments  for  the  fixed  bed  studies  were  carried  out  for 
phenol  and  PCP  at  temperature  from  20°  to  35°C  and  from  35°  to  20°C.   A  fixed 
bed  program  was  modified  by  Mathews  (1984)  to  take  into  account  temperature 
change.   Input  data  of  the  isotherm  parameters,  and  the  kinetic  coeffi- 
cients, Kf  as  well  as  Ds,  were  calculated  from  Equation  4.1  and  from  Equation 
4.6  through  Equation  4.9.   The  time  of  step  change  in  temperature  for  the 
fixed  bed  experiments  was  after  3  hours  and  it  took  30  minutes  to  raise  the 
temperature  from  20°C  to  35°C.   The  input  data  for  the  isotherm  parameters  and 
kf  and  Ds  at  20°C  was  needed  during  initial  195  minutes.   After  that,  the  data 
of  35°C  was  applied  to  predict  breakthrough  curves. 

Figures  4.22  through  4.25  show  the  experimental  and  predicted  data  of  the 
breakthrough  curves  for  phenol  and  PCP.  Here,  also,  the  modified  program  shows 
good  predictions  for  all  cases.   It  is  apparent  that  the  tests  for  step  change 
from  20°  to  35°  show  less  adsorption  than  those  from  35°  to  20°C  from  compa- 
ring Figure  4.22  to  4.23  and  Figure  4.24  to  4.25. 
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Table  4.6 


Operating  Data  for  Fixed  Bed  Experiments  at  20°C 


phenol  PCP 


Total  weight  of  carbon  75.0  gm  75.0  gm 

Column  diameter  3.1  cm  3.1  cm 

Total  bed  height  23.2  cm  23.2  cm 

porosity  of  bed  0.359  0.359 

Influent  concentration  194.12  mg/1  250.20  mg/1 

Total  flowrate  125.00  ml/min  125.00  ml/min 

Surface  loading  0.17  m3/min  m2    0.17  m3/min  m2 

Dg«  number  1188.29  1110.74 

Stb  number  2.98  3.26 

Bic  number  45.11  38.16 


Dg  :  Solute  distribution  parameter  (dimensionless). 
Pa   qe  (l-e)/(e  Co) 

pa   ■  adsorbent  density  including  pore  volume 
e  =  bed  void  fraction 

St  :  Stanton  number  (dimensionless),   kf  t  (l-e)/(R  e  ^) 
T  «  fluid  residence  time  in  fixed  bed 

L/V,  V  :  intersitial  fluid  velocity,  (L/v) 
0  ■  sphericity,  ratio  of  surface  area  of 

equilivalent-volume  sphere  to  actual  surface 

area  of  adsorbent  particle. 


Bi  :  Biot  number  based  on  solid  phase  diffusion 
coefficient  (dimensionless) , 
kf  R  (l-e)/(Dg  Ds  s  0) 
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Table  4.7 

Operating  Data  for  Fixed  Bed  Experiments  at  35°C 

phenol  PCP 

Total  weight  of  carbon  75.0  gut  75.0  gm 

Column  diameter  3.1   cm  3.1   cm 

Total  bed  height  23.2  cm  23.2   cm 

porosity  of  bed  0.359  0.359 

Influent  concentration  195.45  mg/1  243.78  mg/1 

Total  flowrate  125.00  ml/min  125.00  ml/min 

Surface  loading  0.17  m3/min  m2   0.17  m3/min  m2 

Dg  number  918.35  768.39 

St  number  4.06  4.06 

Bi  number  6.73  18.99 
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Table  4.8 

Operating  Data  for  Fixed  Bed  Experiments  : 
Step  Change  from  20°C  to  35°C  after  3  hrs 


phenol 

PCP 

Total  weight  of  carbon 

75.0  gm 

75.0  g» 

Column  diameter 

3.1  cm 

3.1  cm 

Total  bed  height 

23.2  cm 

23.2  cm 

porosity  of  bed 

0.359 

0.359 

Influent  concentration 

194.12  mg/1 

250.20  mg/1 

Total  flowrate 

125.00  ml/min 

125.00  ml/min 

Surface  loading 

0.16  m3/min  m2 

0.16  m3/min  m2 

Dg  number    (20°C) 

1278.42 

976 .  92 

<35°C> 

990.71 

840.58 

St  number     <20°C) 

3.07 

3.35 

(35°C> 

4.27 

5.35 

Bi  number    <20°C) 

29.40 

38.54 

<35°C) 

5.52 

14.83 
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Table  4.9 

Operating  Data  for  Fixed  Bed  Experiments  : 
Step  Change  fro*  35°C  to  20°C  after  3  hrs 

phenol  PCP 

Total  weight  of  carbon  75.0  g»  75.0  ga 

Column  diameter  3.1   cm  3.1   cm 

Total  bed  height  23.2  cm  23.2  cm 

porosity  of  bed  0.359  0.359 

Influent  concentration  195.45  mg/1  243.78  »g/l 

Total  flowrate  125.00  ml/min  125.00  ml/min 
Surface  loading                 0.16  m^/min  m^   0.16  m3/min  m^ 


Dg  number 

(35°C) 

(20°C) 

St  number 

<35°C) 

<20°C> 

Bi  number 

(35°C) 

<20°C> 

999.95 

840.58 

1290.35 

976.92 

4.27 

5.35 

3.07 

3.35 

5.47 

14.83 

39.14 

38.54 
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CHAPTER   5. 


MULTICOMPONENT  ADSORPTION 


5-1  Equilibrium  Studies 

In  the  multiconponent  equilibrium  studies,  experimental  data  of  a  mixture 
of  phenol  and  PCP  were  collected  at  20°C  and  35°C.   Two  models  were  selected 
to  describe  the  bisolute  equilibrium  system.   They  were  : 

1.  Mathews'  model  (1975) 

*i  <Ci/-7i> 

<U  "  -ft. (2.31) 

1  *  SBj  (C/Z-Oj/5-1 

2.  Prausnitz's  model  (IAS  model) 

The  isotherm  parameters.  A,  B,  and  #,  used  in  both  models  are  obtained  from 
the  single  solute  equilibrium  studies. 

In  recent  years,  it  becomes  necessary  to  make  use  of  theoretical  models 
that  allow  the  prediction  of  mixture  equilibria  from  a  minimum  amount  of 
experimental  data.   The  comparison  of  these  two  models  employing  experimental 
data  will  be  discussed  to  find  efficiency  of  their  predictions.   Then,  the 
models  will  be  introduced  to  predict  adsorption  rates  of  multicomponent  system 
which  will  be  discussed  in  the  next  section. 

Tables  5.1  and  5.2  list  the  experimental  and  predicted  values  of  multicom- 
ponent equilibrium  studies  for  Mathews'  model  and  Prausnitz's  model  at  20°C. 
Tables  5.3  and  5.4  show  the  experimental  and  predicted  values  of  both  models 
at  35°C.  The  interaction  coefficients,  171,  are  obtained  from  mixture  equilib- 
rium data  using  single  solute  isotherm  parameters.  A,  B,  and  fi,    developed  for 
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each  solute  species  in  parameter  estimation  procedure.   The  values  of  171  are 
1.015  and  0.4113  for  phenol  at  20°  and  35°C,  respectively;  whereas  for  PCP, 
those  are  0.2095  and  0.1512. 

For  Mathews'  model,  the  percent  deviation  are  26.6,  25.2  for  phenol  and 
21.9,  20.7  for  PCP  at  20°C  and  35°C,  respectively.   For  Prausnitz's  model.  The 
values  are  16.6,  15.8  for  phenol  and  13.1  12.4  for  PCP  at  20°C  and  35°C, 
respectively.   Fro*  The  percent  deviation  standpoint,  it  can  be  seen  that  the 
surface  concentration  predictions  by  Prausnitz's  model  are  in  better  condi- 
tions than  Mathews'  model. 

Both  models  predict  that  PCP  is  adsorbed  more  strongly  than  phenol,  which 
is  evident  from  individul  isotherm.   The  same  results  are  proved  that  adsor- 
ption capacities  are  higher  at  lower  temperature  from  both  models . 

It  is  not  obvious  to  compare  the  results  only  from  the  data  in  Tables  5.1 
to  5.4.   The  plots  of  qexp/qcalc  against  the  solution  concentration  for  both 
models  are  illustrated  in  Figures  5.1  to  5.4  representing  the  Mathews'  and 
Prausnitzs'  models  for  phenol  and  PCP  at  20°  and  35°C.   For  both  temperatures, 
the  prediction  of  Prausnitz's  model  is  seen  to  be  better  than  that  of  Mathews' 
model . 
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Table  5.1 

Multicomponent  Equilibrium  for  Phenol  and  PCP 
at  20°C  by  Mathews's  Equilibrium  Model 


Exp 

Model 

Exp 

Model 

^Phenol 

CPCP 

^Phenol 

^Phenol 

qpcp 

qpcp 

mmoles/1 

noles/qj 

mnol 

es/qm 

4.202 

2.185 

1.5190 

1.4490 

2.6550 

1.9430 

3.277 

1.013 

0.8547 

1.8270 

1.9750 

1.3260 

1.875 

3.843 

E-l 

0.3180 

0.2946 

1.2490 

0.7547 

1.576 

2.914 

E-l 

1.1440 

1.2080 

1.3940 

1 . 1700 

1.218 

1.829 

E-l 

0.9440 

0.9414 

1.2020 

0.7957 

7.550  E-l 

1.386 

E-l 

0.4360 

0.5384 

0.6860 

0 . 5979 

5.190  E-l 

3.945 

E-l 

1 . 1700 

1.0210 

0.4740 

0.3831 

2.680  E-l 

9.799 

E-3 

0.3740 

0.4767 

0.2310 

0.1888 

1.939  E-l 

8.736 

E-3 

0.8800 

0.6602 

1.2350 

0 . 8988 

1.180  E-l 

8.733 

E-3 

0.6100 

0.7290 

0.2387 

0.2861 

1.080  E-l 

5.816 

E-3 

0.7120 

0.6338 

0.5010 

0.4748 

6.315  E-2 

4.983 

E-3 

0.4930 

0.5994 

0.6240 

0.4733 

3.726  E-2 

4.774 

E-3 

0.3090 

0.2734 

0.3330 

0.3395 

2.744  E-2 

4.351 

E-3 

0.4800 

0.3192 

0.1920 

0.2849 

1.898  E-2 

4.253 

E-3 

0.3320 

0.2396 

0.2490 

0.3320 

1.670  E-2 

2.519 

E-3 

0.2490 

0.2247 

0.3570 

0.4287 

1.592  E-2 

1.810 

E-3 

0.2500 

0.1659 

0 . 9900 

0.9717 

3.280  E-3 

9.370 

E-4 

0.0567 

0.5340 

0.6520 

0.5496 
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Table  5.2 

Multicoaponent  Equilibria  for  Phenol  and  PCP 
at  20°C  by  Prausnitz's  Equilibrium  Model 


Exp 

Model 

Exp 

Model 

Cphenol 

CPCP 

"JPhenol 

^Phenol 

qpcp 

qpcp 

] 

■aoles/1 

mmoles/qm 

**oles/qn 

4.202 

2.185 

1.5190 

1.9660 

2.6550 

2.4630 

3.277 

1.013 

0.8547 

0.6549 

1.9750 

1.8700 

1.875 

3.843 

E-l 

0.3180 

0.2934 

1.2490 

1.3807 

1.576 

2.914 

E-l 

1.1440 

1 . 1730 

1.3940 

1.4770 

1.218 

1.829 

E-l 

0.9440 

0.9906 

1.2020 

1 . 3760 

7.550 

E-l 

1.386 

E-l 

0.4360 

0.4560 

0.6860 

0.6083 

5.190 

E-l 

3.945 

E-l 

1 . 1700 

1 . 1380 

0.4740 

0.3376 

2.680 

E-l 

9.799 

E-3 

0.3740 

0.2992 

0.2310 

0.2063 

1.939 

E-l 

8.736 

E-l 

0.8800 

0.9241 

1.2350 

1.3180 

1.180 

E-l 

8.733 

E-3 

0.6100 

0.5152 

0.2387 

0.2087 

1.080 

E-l 

5.816 

E-3 

0.7120 

0.7890 

0.5010 

0.4792 

6.315 

E-2 

4.983 

E-3 

0.4930 

0.4772 

0.6240 

0.5105 

3.726 

E-2 

4.774 

E-3 

0.3090 

0.4236 

0.3330 

0.3220 

2.744 

E-2 

4.351 

E-3 

0.4800 

0.3828 

0.1920 

0.1438 

1.898 

E-2 

4.253 

E-3 

0.3320 

0.3167 

0.2490 

0.2136 

1.670 

E-2 

2.519 

E-3 

0.2490 

0.2938 

0.3570 

0.2919 

1.592 

E-2 

1.810 

E-3 

0.2500 

0.2577 

0.9900 

0.9393 

3.280 

E-3 

9.370 

E-4 

0.0567 

0.0780 

0.6520 

0.4227 
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Table  5.3 

Hulticomponent  Equilibrium  for  Phenol  and  PCP 
at  35°C  by  Mathews's  Equilibrium  Model 


Exp 

Model 

Exp 

Model 

cPhenol 

CPCP 

<JPhenol 

<JPhenol 

qpcp 

qpcp 

■moles/1 

mmoles 

i/gm 

mmol 

es/gm 

3.887 

2.153 

1.5030 

1.0000 

1.6920 

1.4250 

3.009 

1.272 

0.8860 

0.6324 

1.1220 

1.0700 

1.695 

7.363 

E-l 

0.3660 

0.2799 

1.3420 

0.7959 

1.383 

2.972 

E-l 

0.8723 

0.8681 

1.3080 

0.7228 

1.187 

2.714 

E-l 

0.7990 

0.7203 

0.8670 

0.8398 

1.141 

2.443 

E-l 

0.9100 

0.8617 

0.7910 

0.6146 

6.195  E-l 

1.897 

E-l 

0.4510 

0.4783 

1.3350 

0.5742 

4.928  E-l 

1.003 

E-l 

1.0650 

1.0890 

0.8614 

0.7489 

3.265  E-l 

5.163 

E-l 

0.7420 

0.7779 

0.9620 

0.6462 

2.437  E-l 

5.010 

E-2 

0.8290 

0.5934 

0.4260 

0.4889 

1.794  E-l 

3.994 

E-2 

0.5760 

0.7159 

0.7710 

0.6627 

1.791  E-l 

3.876 

E-2 

0.6710 

0.7281 

0.4150 

0.3701 

1.771  E-l 

3.115 

E-2 

0.3740 

0.3300 

0.5500 

0.5133 

9.925  E-2 

2.261 

E-2 

0.2901 

0.3100 

0.5690 

0.2003 

9.963  E-2 

2.210 

E-2 

0.4680 

0.3973 

0.3670 

0.4754 

6.794  E-2 

1.863 

E-2 

0.3140 

0.3464 

0.3890 

0.2979 

4.998  E-2 

1.803 

E-2 

0.2380 

0.2095 

0 . 3400 

0.2920 

2.954  E-2 

1.519 

E-3 

0.4580 

0.2898 

0.5160 

0.5011 

2.442  E-2 

9.675 

E-3 

0.1600 

0.2010 

0.2270 

0.3803 

Table  5.4 

Multiconponent  Equilibrium  for  Phenol  and  PCP 
at  35°C  by  Prausnitz's  Equilibriun  Model 
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Exp 

Model 

Exp 

Model 

Cphenol 

CPCP 

<JPhenol 

^Phenol 

qpcp 

qpcp 

»»oles/l 

ntmoles/gn 

mntol 

es/gm 

3.887 

2.153 

1.5030 

1.3600 

1.6920 

2.1690 

3.009 

1.272 

0.8860 

0.6794 

1 . 1220 

1.3140- 

1.695 

7.363 

E-l 

0.3660 

0.3256 

1 . 3420 

1.2730 

1.383 

2.972 

E-l 

0.8723 

0.8307 

1.3080 

0.9451 

1.187 

2.714 

E-l 

0.7990 

0.8897 

0.8670 

1.3280 

1.141 

2.443 

E-l 

0.9100 

1.2070 

0.7910 

0.5319 

6.195  E-l 

1.897 

E-l 

0.4510 

0.4138 

1.3350 

1.4750 

4.928  E-l 

1.003 

E-l 

1 . 0650 

0.9467 

0.8614 

0.5982 

3.265  E-l 

5.163 

E-l 

0.7420 

0.9301 

0.9620 

1.1950 

2.437  E-l 

5.010 

E-2 

0.8290 

0.7987 

0.4260 

0.3756 

1.794  E-l 

3.994 

E-2 

0.5760 

0.4865 

0.7710 

0.6178 

1.791  E-l 

3.876 

E-2 

0.6710 

0.7756 

0.4150 

0.3167 

1.771  E-l 

3.115 

E-2 

0.3740 

0.3822 

0.5500 

0.4098 

9.925  E-2 

2.261 

E-2 

0.2901 

0.3063 

0.5690 

0.6266 

9.963  E-2 

2.210 

E-2 

0.4680 

0.5307 

0.3670 

0.2093 

6.794  E-2 

1.863 

E-2 

0.3140 

0.4575 

0.3890 

0.2782 

4.998  E-2 

1.803 

E-2 

0.2380 

0.3228 

0.3400 

0.2782 

2.954  E-2 

1.519 

E-3 

0.4580 

0.3002 

0.5160 

0.4092 

2.442  E-2 

9.675 

E-3 

0.1600 

0.2398 

0.2270 

0.1524 
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5-2  Batch  Kinetic  Studies 

Four  runs  with  two  constant  temperatures,  20°C  and  35°C,  and  two  step 
changes,  from  20°C  to  35°C  and  from  35°C  to  20°C,  were  conducted  for  batch 
kinetic  studies  of  multicomponent  solution.   The  operating  conditions  were  the 
same  as  in  single  solute  batch  reactions  as  discussed  in  section  4-2.   Carbon 
size  of  20/25,  stirring  speed  of  700  rpm,  reaction  volume  of  24  liters  solu- 
tion made  up  of  tap  water,  total  carbon  weight  of  6  gm,  and  pH  value  at  7  were 
used  for  all  batch  kinetic  multicomponent  experiments.   The  initial  concentra- 
tion of  the  solutions  included  2.5  x  10~4M  phenol  and  2.5  x  10~4M  PCP. 

Adsorption  rates  for  constant  temperature  were  predicted  by  using  the 
computer  program  developed  by  Mathews  (1975).   The  data  of  both  Mathews'  and 
Prausnitz  's  equilibrium  models  were  used  in  the  computer  program  seperately. 
Figures  5.5  and  5.6  depict  the  results  with  Mathews'  equilibrium  model  at  20°C 
and  35°C;  Figures  5.7  and  5.8  show  the  results  with  Prausnitz's  equilibrium 
model. 

For  Prausnitz's  model  in  predicting  adsorption  rates,  the  percent  devia- 
tion between  the  experimental  and  predicted  data  are  2.58  for  phenol  and  16.89 
for  PCP  at  20°C  as  well  as  5.82  for  phenol  and  19.85  for  PCP  at  35°C.   For 
Mathews'  model,  the  values  are  4.95  for  phenol  and  9.43  for  PCP  at  20°C  as 
well  as  6.00  for  phenol  and  10.62  for  PCP  at  35°C.   Both  models  gave  reason- 
able predictions.    It  seems  that  adsorption  rate  model  using  Prausnits's 
equilibrium  correlation  does  not  give  as  good  a  prediction  as  the  Mathews' 
model. 

For  step  change  experiments,  adsorption  rates  were  predicted  by  using  rate 
model  incorporating  the  Prausnitz's  equilibrium  models.   A  modification  of  the 


97 


computer  program  incorporating  Prausnitz's  equilibrium  model  developed  by 
Mathews  (1984)  was  used  to  predict  kinetic  curves  for  change  in  temperature, 
similar  to  previous  sections.   The  result  of  adsorption  rate  for  step  change 
from  20°C  to  35°C  is  represented  in  Figure  5.9.   The  percent  deviations  bet- 
ween the  predicted  and  experimental  data  are  4.29  for  phneol  and  20.32  for 
PCP.   It  took  20  minutes  for  the  temperature  of  solution  to  drop  from  35°C  to 
20°C.   In  predicting  adsorption  rates  from  35°C  to  20°C,  the  result  give  as 
reasonable  agreement  as  can  be  seen  from  Figure  5.10.   The  percent  deviations 
are  9.90  for  phenol  and  9.18  for  PCP.   The  possible  explanation  of  the 
difference  between  experimental  and  predicted  data  is  that  the  procedure  of 
adding  dry  ice  to  decrease  the  temperature  affect  the  stirring  speeds  and 
induces  vaporization  of  adsorbates. 
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5-3  Fixed  Bed  Studies 

Two  runs  were  made  in  the  fixed  bed  experiments  with  multicomponents 
solutions,  i.e.,  step  changes  from  20°C  to  35°C  and  from   35°C  to  20°C  after  3 
hours.   The  total  flowrate  and  total  weight  of  carbon  used  in  two  runs  were 
the  same  as  those  in  the  single  solute  studies.   The  breakthrough  curves  of 
two  runs  were  obtained  from  conducting  multicomponent  solution  consisting  1  x 
10_3M  of  phenol  and  2  x  10"3M  of  PCP  as  initial  concentration.   The  runs  were 
terminated  till  the  C/Co  values  of  PCP  approaching  0.95.   The  duration  of 
experiments  was  12  hours. 

Figure  5.11  shows  the  experimental  data  of  the  step  change  experiments 
from  20°  to  35°C;  while  Figure  5.12  depicts  the  data  from  35°  to  20°C.   The 
reason  for  C/C0  values  of  phenol  greater  than  1  is  because  of  displacement  of 
previously  adsorbed  phenol  by  PCP  and  elution  of  phenol  into  liquid  phase. 
The  amount  of  phenol  displaced  by  PCP  reached  maximum  at  the  highest  points 
shown  in  the  plots.   After  the  peak,  the  displacement  became  slower  and  the 
values  of  C/Co  of  phenol  decreasd  as  time  progressed. 

From  the  plots,  it  can  be  seen  that  for  phenol,  the  C/Co  values  for  step 
change  from  35°C  to  20°C  were  greater  than  those  from  20°C  to  35°C;  while  for 
PCP,  the  values  from  35°C  to  20°C  were  less  than  those  from  20°C  to  35°C.   The 
mechanisms  and  phenomenon  of  the  multicomponent  fixed  bed  adsorption  for  step 
change  in  temperature  are  very  complex  because  the  time  for  changing  tempera- 
ture is  also  a  variable  affecting  adsorption  profiles.   The  observation  of 
variables  on  multicomponent  fixed  bed  adsorption  is  not  within  the  scope  of 
this  research.   Further  investigations  are  recommended  on  this  Doint, 
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CHAPTER   6. 
CONCLUSIONS  AND  RECOMMENDATIONS 

6-1  Conclusions 

1.  Temperature  affects  both  adsorption  capacity  fro*  equilibrium  experiments 
and  adsorption  rates  from  batch  kinetic  experiments. 

2.  The  three  parameter  equation  described  adsorption  capacity   well  for  both 
phenol  and  PCP  at  four  different  temperatures.   The  extended  three  parame- 
ter equation  vas  developed  to  show  the  temperature  effect  on  capacity 
based  on  data  at  four  temperatures. 

3.  The  homogeneous  solid  phase  diffusion  model  described  satisfactorily 
adsorption  rates  of  single  solute  with  constant  temperature.   The  values 
of  kf  and  Ds  obtained  from  the  batch  kinetic  model  were  correlated  with 
temperature  based  on  the  Arrhenius'  law.   The  modified  kinetic  model  was 
tested  for  runs  with  step  changes  from  20°C  to  35°C  and  35°C  to  20°C  and 
showed  good  copoblity  for  predicting  adsorption  rate. 

4.  The  fixed-bed  program  based  on  the  homogeneous  solid  phase  diffusion  model 
was  used  to  predict  breakthrough  curves  of  phenol  and  PCP  at  20°C  and 
35°C.   The  operation  of  thea  fixed  bed  at  20°C  was  found  to  be  better  than 
at  35°C  from  the  breakthrough  curves.   The  modified  fixed-bed  program  used 
in  predicting  breakthrough  curves  for  step  changes  was  found  to  provide 
excellent  prediction. 
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5.  Both  Mathews'  model  and  Prauanitz'a  model  gave  reasonable  estimation  for 
multicomponent  equilibria.   However,  accurcy  of  multicomponent  rate  model 
predictions  seems  to  be  limited  by  the  inability  to  predict  mixture  equi- 
libria satisfactorily. 

6.  In  the  adsorption  rates  studies  for  multicomponent  solution,  the  adsor- 
ption rate  model  based  on  IAS  model  provide  good  prediction  for  runs  with 
constant  temperature.   In  addition,  the  rate  model  which  was  modified  to 
include  temperature  changes  gave  predictions  in  moderately  good  agreement 
for  runs  with  step  change  in  temperature. 

7.  Experimental  results  from  multicomponent  breakthrough  curves  showed  that 
the  C/Co  values  of  phenol  for  step  change  from  35°C  to  20°C  were  greater 
than  those  of  from  20°C  to  35°C;  while  the  values  of  PCP  for  step  change 
from  35°C  to  20°C  were  less  than  those  of  from  20°C  to  35°C. 
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6-2  Recommendations 

1.  Extend  the  application  of  the  Modified  three  parameter  equation,  which  is 
available  for  several  temperature,  to  other  priority  pollutants  and  test 
with  other  ranges  of  temper taures . 

2.  The  rate  model  and  the  fixed-bed  model  were  found  in  good  agreement  in 
predicting  adsorption  rates  and  breakthrough  for  a  single  solute  of  phneol 
and  PCP.   Extend  the  application  of  the  models  and  test  the  efficiency 
with  other  priority  pollutants. 

3.  The  modified  models  for  studing  step  change  can  be  extended  to  take  into 
account  time  as  a  variable,  i.e.,  changing  temperature  at  different  times 
to  test  the  efficiencies  of  the  models  in  predicting  batch  kinetic  data  or 
breakthrough  curves.    The  the  model  can  be  verified  by  running  experi- 
ments and  calculating  predictive  data,  which  were  conducted  with  testing 
different  priority  pollutants,  changing  temperature  at  different  time,  and 
step  change  with  different  temperature. 
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ABSTRACT 

A  mathematical  model  of  fixed  bed  adsorption  was  modified  to  take  into 
account  temperature  as  a  variable.  The  model  was  employed  to  compare  the 
performance  of  fixed  beds  at  constant  temperature  and  step  changes  in  tempera- 
ture. The  equations  which  include  the  kinetic  coefficients,  film  transfer  and 
solid  phase  diffusion  in  particle,  and  adsoprtion  capacity  as  variables  were 
solved  by  orothogonal  collocation  methods. 

Phenol  and  parachlorophenol  (PCP)  were  used  as  adsorbates.   Four  constant 
temperatures,  20°,  25°,  30°.  and  35°C  .  as  well  as  two  temperature  step 
change,  from  20°  to  35°C  and,  from  35°  to  20°C  were  tested.   For  single 
solute,  isotherm  at  four  temperatures  were  conducted.  An  equation  is  developed 
to  show  temperature  effect  on  capacity  based  on  data  at  four  different  tempe- 
ratures.  Batch  kinetic  studies  were  carried  out  at  four  temperatures  and  two 
step  changes.   Film  transfer  and  solid  phase  diffusion  coefficients  obtained 
from  batch  kinetic  model  were  correlated  with  temperature.   The  correlated 
film  transfer  and  solid  phase  diffusion  coefficients  were  used  to  predict 
adsorption  rates  for  step  changes  in  temperature.  For  multicomponent,  equilib- 
rium and  batch  kinetic  studies  were  accomplished  for  temperatures  of  20°  and 
35°C.   Hathews'  and  Prausnitz's  models  were  adopted  to  fit  the  experimental 
data.   But  only  Prausitz's  model  was  modified  to  predict  the  batch  runs  with 
step  change  because  of  its  more  accuracy  in  describing  constant  temperature 
tests. 

In  fixed  bed  runs,  experimental  data  were  tested  against  the  mathematical 
model  with  20°  and  35°C.   For  step  change  tests,  the  modified  model  using 


variables  from  extended  isotherm  and  correlated  Kf,  D3  was  applied  in  predic- 
ting breakthrough  profile.  Multiconponent  adsorption  in  fixed  beds  for  step 
change  in  temperature  were  conducted  to  compare  breaktnrough  profile  of  each 
solute. 


